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Microbial  fuel  cell  (MFC)  is  a  new  approach  to  accomplish  simultaneous  waste  treatment  and  alter¬ 
native  energy  production.  An  MFC  employs  exoelectrogenic  bacteria  to  directly  convert  biodegradable 
substrates  to  electricity.  The  electricity  production  from  MFCs  can  be  maximized  at  a  certain  operating 
condition,  which  is  defined  as  the  maximum  power  point  (MPP).  However,  it  is  difficult  to  maintain 
an  MFC  at  its  MPP  during  operation,  because  the  power  output  from  MFCs  fluctuates  constantly  due  to 
the  changes  of  microbial  activities  or  environmental  conditions.  This  study  developed  a  maximum  power 
point  tracking  (MPPT)  technique  using  digitally  controlled  potentiometers  for  hysteresis  controller  based 
energy  harvesting  system.  Results  with  lab  scale  MFC  reactors  have  shown  that  the  system  can  track  the 
MPP  and  maintain  the  maximum  energy  harvesting  in  real-time,  suggesting  that  the  proposed  system 
can  capture  the  most  available  energy  from  MFCs  at  the  most  efficient  operating  point. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Microbial  fuel  cell  technology  has  been  intensively  researched 
in  recent  years  as  a  novel  technology,  because  it  offers  a  solu¬ 
tion  for  energy  and  environmental  sustainability  by  simultaneously 
performing  waste  treatment  and  generating  electricity  -  a  clean, 
distributed,  and  even  remote  energy  source.  MFCs  use  electrochem- 
ically  active  bacteria  to  extract  the  electrons  from  biodegradable 
substrates  in  the  anode  chamber  and  transfer  them  to  the  anode 
electrode.  These  electrons  flow  from  the  anode  through  an  exter¬ 
nal  circuit  to  the  cathode,  where  they  reduce  an  electron  acceptor 
such  as  oxygen  or  ferricyanide.  The  reported  maximum  power 
density  from  lab  scale  air-cathode  reactors  has  been  increased  by 
orders  of  magnitude  in  the  last  few  years,  from  less  than  1  mW  m-2 
to  6.9  Wm-2  [1,2].  The  typical  MFC  power  density  measurement 
method  involves  operating  the  reactor  with  a  static  external 
resistance  or  applied  potential,  and  then  transiently  obtaining 
polarization  data  by  applying  a  series  of  different  external  resis¬ 
tances  or  performing  an  electrochemical  scan  using  a  potentiostat 
[3,4].  The  maximum  power  density  is  achieved  when  the  applied 
external  resistance  is  equal  to  the  MFC  internal  resistance  [3-5]. 
Such  characterizations  represent  the  potential  of  MFC  power  out¬ 
put,  but  they  don’t  measure  the  actually  usable  power  output, 
because  the  generated  electricity  is  dissipated  into  heat  instead  of 
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being  utilized  by  electronics.  Moreover,  the  fixed  resistance  can¬ 
not  recover  the  maximum  power  output  during  MFC  operation, 
because  the  internal  resistance  of  an  MFC  varies  with  changes  in 
operational  parameters,  such  as  substrate  concentration,  pH,  and 
temperature  [6-8]. 

For  efficient  harvesting  and  usage  of  the  MFC  energy,  tech¬ 
nologies  capable  of  real-time  maximum  power  point  tracking 
(MPPT)  and  the  maximum  available  power  capturing  need  to  be 
developed.  MFCs  may  lost  more  than  50%  of  produced  power 
across  the  internal  resistance  if  the  operating  voltage  is  not  at  the 
maximum  power  point  (MPP)  voltage  [9].  In  addition,  compared 
to  other  physiochemical-reaction  based  processes,  MFC  perfor¬ 
mance  depends  on  microbial  activity,  which  is  relatively  slow 
and  susceptible  to  many  environmental  condition  changes  such  as 
instantaneous  output  power  level,  accumulated  extracted  energy, 
bacteria  community  and  activity  shifts.  Furthermore,  the  output 
voltage  of  MFC  does  not  remain  constant  as  the  current  output 
increases,  making  the  real-time  energy  tracking  and  capturing 
more  challenging. 

Popular  MPPT  techniques  such  as  perturbation  and  observation 
(P&O)  or  gradient  method  for  photovoltaic  systems  and  hydrogen 
fuel  cells  were  introduced  to  MFC  systems  for  real-time  optimal 
selection  of  the  external  resistance  [5,9].  These  techniques  have 
been  utilized  to  identify  the  value  of  proper  external  resistance  to 
operate  the  MFC  at  MPP.  However,  finding  the  optimal  external 
resistance  still  cannot  capture  and  utilize  the  energy  generated  by 
MFCs.  Another  approach  for  harvesting  system  is  using  capacitor- 
based  circuits  such  as  super  capacitors  and  charge  pumps,  which 
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Energy  Harvesting  Converter  Voltage  Boost  Converter 


Fig.  1.  Block  diagram  of  hysteresis  controller  based  MFC  energy  harvesting  system. 


capture  the  energy  passively  and  transfer  it  to  a  boost  converter 
[10-12].  But  these  systems  are  not  suitable  for  MPPT  either  because 
of  the  limited  controllability  on  operating  condition. 

In  this  study,  an  integrated  control  system  that  can  perform  both 
real-time  MPPT  and  the  maximum  available  power  harvesting  has 
been  developed.  The  system  does  not  utilize  charge  pumps  and  was 
operated  under  batch  conditions  with  a  hysteresis  controller  based 
energy  harvesting  system.  The  power  output  tracked  and  regulated 
by  the  controller  was  compared  with  the  power  generation  from 
fixed  external  resistance  operation,  and  the  real-time  parameter 
changing  was  compared  with  the  polarization  power  curve. 

2.  Materials  and  methods 

2.1.  Hysteresis  controller  based  MFC  energy  harvesting  system 

A  hysteresis  circuit  based  MFC  energy  harvesting  system  has 
been  utilized  in  this  paper.  The  energy  harvester  is  the  front-end  of 
two  layers  of  DC/DC  converters:  the  first  converter  (energy  har¬ 
vester)  captures  the  energy  from  an  MFC,  and  the  second  layer 
converter  (voltage  booster)  increases  the  voltage  to  a  suitable  level 
for  the  load.  Fig.  1  shows  the  block  diagram  of  the  overall  system, 
and  operation  details  can  be  found  in  authors’  previous  work  [13]. 

The  system  maintains  the  operating  voltage  at  the  MPP  of 
MFC  by  switching  the  harvesting  converter  switches  Sn  and  Si2. 
The  operation  of  the  energy  harvester  in  the  first  layer  consists 


of  two  modes,  CFIARGE  and  DISCHARGE:  during  CHARGE  mode, 
the  energy  harvester’s  main  switch  Sn  is  turned  on  to  get  the 
energy  extracted  from  the  MFC  and  stored  in  the  inductor  L\. 
The  MFC  voltage  VMfc  decreases  due  to  the  increasing  output 
current.  This  switch  is  off  during  the  DISCHARGE  mode,  and  the 
energy  in  the  inductor  L\  is  discharged  to  the  storage  capacitor  C\ 
through  the  secondary  switch  Si2.  The  MFC  voltage  recovers  as  cur¬ 
rent  decreases.  The  switches  are  toggled  complimentarily  when 
the  MFC  voltage  reaches  lower  threshold  and  upper  threshold  in 
the  CHARGE  and  DISCHARGE  mode,  respectively.  Generally,  an  N- 
channel  MOSFET  is  utilized  for  Sn  because  of  high  switching  speed 
and  low  conduction  loss.  For  Si2,  a  diode  can  be  used  for  simple 
circuit  configuration,  but  it  suffers  from  high  forward  voltage  drop. 
For  better  efficiency,  a  synchronous  converter  can  be  implemented 
using  a  P-channel  MOSFET  as  Si2,  but  it  requires  a  complex  gate 
drive  circuitry  because  it  is  a  floating  switch.  The  upper  and  lower 
voltage  thresholds  are  given  as  VthH  =  VCc  -^2/(^2  +  (^i//^3)}  and 
VthL  =  VCc  •  (R2//R3 )/ {Ri  +  (R2//^3  )}•  The  frequency  of  power  extrac¬ 
tion,  i.e.,  the  transition  between  CHARGE  and  DISCHARGE  mode, 
is  automatically  controlled  according  to  the  MFC’s  condition  so 
that  the  MFC  voltage  can  be  confined  in  a  pre-defined  voltage 
band,  VthH  -  VthL,  and  at  the  same  time  enough  recovery  time  of 
the  MFC  reactor  is  ensured.  This  scheme  is  especially  effective  for 
the  weak  MFC  units  which  experience  significant  internal  volt¬ 
age  drop  with  an  output  current  increase.  The  second  layer  with 
a  standard  DC/DC  boost  converter  increases  the  output  voltage  to 


Fig.  2.  Proposed  MPPT  controller  using  digital  potentiometers. 
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a  practically  usable  level,  but  it  can  be  optional  if  the  load  takes 
burst-type  power. 


2.2.  Hysteresis  controller  based  maximum  power  tracking  and 
harvesting 


In  order  to  track  and  maintain  the  operating  point  at  MPP 
continuously  despite  the  conditions  (and  therefore  polarization 
curve)  changing  during  MFC  operation,  digitally  controllable  poten¬ 
tiometers  have  been  integrated  into  the  hysteresis-based  energy 
harvester  (Fig.  2).  The  MFC  voltage  and  current  outputs  are  mon¬ 
itored  and  fed  into  the  MPPT  controller  for  instantaneous  output 
power  calculation.  Then  the  MPPT  controller  generates  control  sig¬ 
nals  for  R2  and  R3  to  adjust  the  operating  voltage  based  on  the 
calculated  power.  A  reduced  degree-of-freedom  control  scheme 
was  created  for  the  hysteresis  controller  to  overcome  the  coupling 
problem  of  the  multiple  resisters.  The  system  controls  the  volt¬ 
age  thresholds  and  hysteresis  band  based  only  on  R2.  The  feedback 
resistor  R3  is  adjusted  as  a  function  of  R2. 

The  upper  and  lower  voltage  thresholds  can  be  written  as 


VthH  =  VCC  • 


^1^2  +  ^2^3 
i?2  +  ^2^3  +  R]R3 


VthL  =  Vcc  • 


R2R3 

i?2  +  ^2^3  +  P1P3 


(1) 

(2) 


If  R3  is  adjusted  as  a  dependent  variable  of  R2,  such  as  R3=I<  R2, 
the  thresholds  can  be  given  as  follows. 

V  -V  R1  +  ^R2  fos 

ChH  CC  (/<■  _|_  1  )R1  +  KR2  ^ 

VthL  =  Vcc  ■  (K  +  1  )R,  +  KR2  **  Vcc  '  R,  +  R2  ^ 

Then  the  hysteresis  voltage  band  becomes 

Vband  =  Vcc  '  +  1  )Rl  +  KRl  %  ^CC  K  R]  +  R2  '  ^ 

It  has  been  shown  in  (3)  that  the  lower  threshold  can  now 

be  controlled  with  fixed  R 1  and  variable  R2,  and  the  coefficient 
I<  determines  the  hysteresis  voltage  band  with  given  R 1  and  R2 
in  (5).  Although  the  hysteresis  voltage  band  varies  as  the  voltage 
thresholds  change,  it  can  be  maintained  reasonably  small  by  using  a 
simple  linear  coefficient,  such  as  I<=  C{1  -  R2/R^ },  where  C  is  a  con¬ 
stant.  For  example,  the  hysteresis  voltage  band  will  be  18-20  mV  if 
C  =  50  and  1 .5  V  reference  voltage. 

For  the  experiment  in  this  paper,  the  P&O  method  was  imple¬ 
mented.  A  discrete-domain  equation  for  R2  adjustment  based  on 
P&O  technique  can  be  given  as  follows,  where  P0  is  the  calculated 
output  power  and  I<R  is  the  adjustment  gain  to  control  the  step  size. 


R2[n  +  ]]  =  R2[n]  +  I(R 


P0[n]  -  P0[n  -  1] 
R2[n]-R2[n-  1] 


(6) 


The  controller  adjusts  R2  in  the  direction  that  increases  the  out¬ 
put  power  and  finds  the  instantaneous  MPP. 


3.  Experimental 

3.1.  MPPT  controller 


For  MPP  tracking,  the  MPPT  controller  block  in  Fig.  2  has 
been  implemented  using  digitally  controllable  potentiometers. 
Microchip  MCP4251  5k£2  and  100  kQ  has  been  chosen  for  R2  and 
P3,  respectively.  Each  MCP4251  has  two  potentiometer  channels 
and  for  P3,  they  are  connected  in  series  to  double  the  resistance 
range.  MCP4251s  are  controlled  using  serial  peripheral  interface 
(SPI)  bus  from  TI  320F28335  Digital  Signal  Processor  (DSP)  control 
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Fig.  3.  MPPT  experiment  setup. 


board.  The  voltage  and  current  measurement  and  power  calculation 
routines,  and  resistance  control  algorithm  have  been  implemented 
in  DSP  firmware.  In  this  experiment,  R2  and  R3  were  adjusted  as 
follows: 

A R2  =  sgn(AP0  x  AR2o)^rjT— —  (7) 

R2[n  +  l]  =  R2[n]  +  AR2  (8) 

R3[n  +  1]  =  50 ■  (i_R2[”  +  1])  - R2 [ n  +  1  ]  (9) 

where,  AP0  =  P0[n] -P0[n -  1],  AR2o  =  R2[n]-R2[n-\]  and  sgn(x) 
function  returns  the  polarity  of  x.  A  fixed  resistance  step  and  unity 
adjustment  gain  was  used  in  this  experiment.  The  Nstep  is  the  max¬ 
imum  number  of  the  digital  potentiometer  steps.  MCP4251  has 
257  steps,  which  makes  the  R2  step  size  19.45  £2.  The  R2  has  been 
adjusted  one  step  size  up  or  down  according  to  the  polarity  of 
AP0  x  AR2o  at  every  control  interval. 

To  calculate  the  instantaneous  power  P0,  the  MFC  output  cur¬ 
rent  has  been  measured  using  Tektronics  A622  current  probe  and 
fed  into  the  DSP  A/D  converter  as  well  as  the  MFC  voltage.  Due  to  the 
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very  low  output  current  level,  a  100-turn  coil  and  a  non-inverting 
amplifier  with  x21  gain  has  been  used  for  current  measurement. 
The  output  voltage  boost  converter  in  second  layer  of  Fig.  1  was 
omitted  in  this  MPPT  experiment.  The  MPPT  capable  energy  har¬ 
vesting  converter  and  the  MFC  reactor  used  in  the  experiment  are 
shown  in  Fig.  3. 

3.2.  MFC  reactor  and  operation 

Two-chamber  cubic-shaped  MFC  reactor  was  constructed  and 
the  diameter  of  each  chamber  was  7  cm  [  1 1  ].  The  anode  and  cathode 
chambers  were  separated  by  a  cation  exchange  membrane  (CEM, 
CMI  7000,  Membranes  international,  NJ)  and  the  volume  of  each 
chamber  was  140mL.  Fleat  treated  graphite  brushes  were  used 
as  the  anodes,  and  plain  carbon  paper  was  selected  as  the  cath¬ 
ode  material  [14].  Phosphate  (50  mM)  buffered  ferricyanide  was 
used  at  the  catholyte  to  minimize  the  cathode  effects  on  system 
performance.  After  inserting  the  electrodes,  the  volume  of  each 
chamber  was  approximately  140mL.  Reactors  were  inoculated 
using  effluent  from  another  MFC  operated  in  the  same  lab  for  more 
than  one  year.  Anolyte  medium  solution  was  prepared  containing 
1.25 gL_1  CH3COONa,  03\gl~A  NH4C1,  0A3gL~'  KC1,  3.321  gL”1 
NaH2PO4-2H2O,10.317gL-1  Na2HP0412FI20, 12.5  mLL-1  mineral 
solution,  and  5  mLL-1  vitamin  [15,16].  Reactors  were  operated 
in  fed-batch  mode  at  room  temperature  and  refilled  with  new 
medium  solution  as  required  by  electrical  characterization.  Polar¬ 
ization  test  was  performed  by  a  potentiostat  (Gamry  Instruments, 
PA,  USA)  at  a  scan  rate  of  0.1  mV  s-1 ,  and  the  power  curve  and  MPP 
can  be  seen  in  Fig.  4. 

4.  Results  and  discussion 

Fig.  5  clearly  shows  the  real-time  operating  point  control  and 
tracking  of  the  MPP  on  the  polarization  curve  of  MFC  under  test. 
The  controlled  operating  voltage  and  power  density  (bold)  follows 
the  polarization  curve  until  it  reaches  the  MPP.  It  can  be  seen  in 
Fig.  6  that  the  proposed  system  tracks  the  MPP  that  is  changing 
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Fig.  5.  MPPT  control  experiment  result  on  polarization  curve. 


Fig.  6.  MPPT  control  experimental  result  on  time  axis. 


Fig.  7.  Detailed  operation  of  energy  harvesting  controller.  Top:  on/off  status  of  switch  Sn.  Middle:  MFC  voltage.  Bottom:  MFC  current. 
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Fig.  8.  Width  of  the  hysteresis  voltage  band  of  constant  coefficient  I<=  50  and  linear 
coefficient  K=  50{1  -  0.5(R2/Ri  )}. 


over  time,  which  has  the  result  of  MPPT  control  experiment  plot¬ 
ted  on  time  axis.  The  MPPT  algorithm  was  activated  when  the  MFC 
voltage  is  0.52  V  at  time  zero  and  adjusted  R2  at  every  10  s  until 
the  output  power  changes  its  direction.  Then  the  power  calcula¬ 
tion  and  resistance  adjustment  interval  was  increased  to  60  s.  The 
controller  found  the  MPP  at  around  0.32  V  at  around  460  s,  and 
maintained  the  operating  point  at  the  changing  MPP  afterwards 
by  adjusting  operating  voltage.  It  can  also  be  seen  that  the  output 
current  is  maintained  at  a  constant  level  but  the  voltage  and  power 
droop  as  the  MFC  generates  power  continuously.  An  example  of 
detailed  operation  of  the  harvesting  controller  is  shown  in  Fig.  7. 
The  top  waveform  shows  the  on/off  status  of  switch  Sn.  The  sys¬ 
tem  is  in  CHARGE  and  DISCHARGE  mode  when  the  switch  is  on  and 
off,  respectively.  It  can  be  seen  that  the  MFC  voltage  is  increased 
and  decreased  in  each  mode  in  second  waveform.  At  this  particular 
instance,  the  MFC  output  voltage  is  controlled  at  around  320  mV. 
The  MFC  current  is  shown  in  the  last  waveform. 

Unlike  many  abiotic  power  systems  whose  power  output  is 
considered  constant  due  to  stable  physical  or  chemical  reac¬ 
tions,  the  output  power  of  the  MFC  reactor  is  determined  by  the 
dynamic  microbial  activity  in  the  anode  chamber.  The  bacterial 
community  and  activity  vary  constantly  due  to  the  changes  of 
microbial  metabolism  and  environmental  conditions  such  as  pH, 
temperature,  or  substrate  concentration,  so  the  associated  power 
production  changes  accordingly.  Such  variations  can  be  minimized 
by  optimizing  microbial  growth  conditions  through  system  oper¬ 
ation  to  stabilize  environmental  factors,  but  the  stabilization  of 
microbial  activity  generally  takes  hours  or  longer.  Such  character¬ 
istics  are  reflected  in  the  voltage  waveform  in  Fig.  6,  where  the 
voltage  declined  gradually  due  to  the  slow  response  of  microbial 
activity.  Substrate  depletion  in  batch  operation  may  also  cause 
gradual  voltage  decline,  and  continuous  flow  operation  will  be  fur¬ 
ther  conducted  to  avoid  such  effects.  To  our  best  knowledge,  no 
control  system  has  been  reported  that  is  able  to  track  the  MPP 
in  situ  while  capturing  actually  usable  energy  from  such  microbial 
systems. 

As  can  be  seen  in  (5),  the  hysteresis  voltage  band  is  a  func¬ 
tion  of  resistances  R\ ,  R2  and  coefficient  I<.  Although  the  control 
is  simplified  by  decreasing  the  degree-of-freedom,  the  width  of  the 
hysteresis  band  varies  as  the  ratio  of  R2  to  R !  changes.  If  a  constant 
value  C  is  used  for  I<  over  the  zero  to  one  range  of  R2  to  Ri  ratio, 
the  initial  voltage  band  when  R2  =  0  is  given  as  1/C  V  and  the  maxi¬ 
mum  voltage  band  deviation  is  almost  50%  when  R2/R\  =  1 .  This  can 


Fig.  9.  Operating  voltage  thresholds  according  to  R2  variation.  The  horizon¬ 
tal  axis  denotes  relative  R2  resistance  to  R\.  A  1.5  V  voltage  reference  and 
JO  50{1  -  0.5(K2/Ki  )}  is  used. 


result  in  inconsistent  operational  conditions.  This  band  width  varia¬ 
tion  can  be  reduced  with  variable  coefficient  I<.  If  I< changes  linearly 
as  a  function  of  the  resistance  ratio,  such  as  I<=C{  1  -0.5(R2/Ri)}, 
the  maximum  variation  of  the  band  width  is  reduced  to  about  11%, 
so  that  the  reasonably  constant  operating  condition  can  be  main¬ 
tained  over  a  wide  voltage  range.  In  Fig.  8,  the  hysteresis  voltage 
band  by  the  constant  coefficient  I<=C=  50  and  linear  coefficient 
K  =  50{1  -0.5(R2/Ri)}  are  shown.  It  can  be  seen  that  the  hystere¬ 
sis  voltage  band  variation  with  linear  coefficient  I<  is  significantly 
smaller  compared  to  the  one  with  the  constant  coefficient  K.  Fig.  9 
shows  the  upper  and  lower  voltage  threshold  with  linear  coefficient 
/C=50{1-0.5(R2/Ri)}. 

A  general  DSP  microcontroller  board  has  been  utilized  for  this 
experiment  and  the  power  for  the  control  circuit  has  been  supplied 
from  an  external  power  supply.  It  is  because  the  controller  is  not 
designed  for  the  low-power  application  and  the  minimal  power 
consuming  control  circuit  design  has  not  been  investigated  in  this 
paper.  However,  using  low  power  parts  and  a  CMOS  ASIC  would  be 
a  feasible  approach  for  minimal  power  consuming  control  circuit 
implementation  and  self-sustainability. 

5.  Conclusions 

In  this  paper,  a  hysteresis  controller  based  MFC  energy  harvest¬ 
ing  system  that  is  capable  of  maximum  power  point  tracking  and 
maximum  available  energy  harvesting  has  been  developed  using 
digitally  controlled  potentiometers  for  the  first  time.  Compared  to 
the  traditional  systems,  the  proposed  system  eliminated  external 
resistances  and  enabled  simultaneous  MPP  tracking  and  maximum 
MFC  energy  capture  in  real-time.  This  approach  provides  a  new  way 
of  robust  and  efficient  MFC  operation  for  maximizing  output  power 
and  carries  potential  to  significantly  advance  the  development  and 
application  of  microbial  fuel  cell  technology. 
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